

/*

Table B18

*/



local nametable = "tables/rev_probgrowinglowhighBorderFE.tex"




*-----
eststo clear
local column_number = 1
local column_name = "column" + "`column_number'"
*---



	
foreach varg in "lowg" "higg"{
	forvalues xd=1(1)2{
	
		if `xd'==1{
			local xdist = 50
		}
		if `xd'==2{
			local xdist = 25
		}



		use "public_data/g30.dta", clear
		
replace groupfe = border

		drop if adist>`xdist'		

		
		
		gen hign 	= 0
		gen lown 	= 1
		gen higc = " "
		gen lowc = " "
				
		foreach crop_areat of varlist *_areat{
		
			capture drop inter*

			*defines a variable that measures whether crop is cultivated
			qui bys groupfe: egen interarea = max(`crop_areat')
			
			local crop 	= subinstr("`crop_areat'","_areat","",.)
			if "`crop'"=="cassavayam"{
				local crop = "cassava"
			}
			if "`crop'"== "cocoacoffeetea"{
				local crop = "cocoa"
			}
			if "`crop'"== "wspotato"{
				local crop = "spotato"
			}
			if "`crop'"== "bancocon"{
			local crop = "coconut"
			}
			if "`crop'"== "vegetables"{
				local crop = "tomato"
			}
			if "`crop'"== "pulses"{
				local crop = "phaseolus"
			}
			if "`crop'"== "cerealoth"{
				local crop = "barley"
			}	
			local crop_null 				= "`crop'" + "_null"
			
			
*			bys groupfe: egen intermed = mode(`crop_null'), minmode			
gen intermed = `crop_null'
			qui replace hign  = intermed 		if intermed>=hign&interarea>0
			qui replace lown  = intermed 		if intermed<=lown&interarea>0
			qui replace higc  = "`crop'" 		if intermed>=hign&interarea>0
			qui replace lowc  = "`crop'" 		if intermed<=lown&interarea>0
			capture drop inter*
		}

		*Computes probabilty that highest or lowest risk failure crop cultivated
		gen higg = 0
		gen lowg = 0
		
		foreach crop_areat of varlist *_areat{
		local crop 		= subinstr("`crop_areat'","_areat","",.)
		if "`crop'"=="cassavayam"{
			local crop = "cassava"
		}
		if "`crop'"== "cocoacoffeetea"{
			local crop = "cocoa"
		}
		if "`crop'"== "wspotato"{
			local crop = "spotato"
		}
		if "`crop'"== "bancocon"{
		local crop = "coconut"
		}
		if "`crop'"== "vegetables"{
			local crop = "tomato"
		}
		if "`crop'"== "pulses"{
			local crop = "phaseolus"
		}
		if "`crop'"== "cerealoth"{
			local crop = "barley"
		}
		local crop_null = "`crop'"+"_null"
		qui replace higg  = 1 if `crop_null'==hign&`crop_areat'>0
		qui replace lowg  = 1 if `crop_null'==lown&`crop_areat'>0
		}

		gen cropn 	  = higc if "`varg'"=="higg"
		replace cropn = lowc if "`varg'"=="lowg"

		replace weight = 1
		
		gen interdep = `varg'


		*computes mean of dependent var
		egen interm = mean(interdep)
		local vmeandep = interm[1]
		local meandep = string(`vmeandep')
		if index("`meandep'",".")==1{
			local meandep = "0"+substr("`meandep'",1,3)
		}
		if index("`meandep'",".")<5&index("`meandep'",".")>1{
			local meandep = substr("`meandep'",1,4)
		}
	

		*computes the number of country pairs
		egen intercp = group(border)
		egen intercpmax = max(intercp)
		local countrypair = intercpmax[1]
		drop intercp*			
	
	
*		xi: areg interdep   varofi   [aw=weight],a(groupfe) cluster(border)
		reghdfe interdep varofi [aw=weight], absorb(groupfe) vce(cluster groupfe country)
		eststo `column_name'
		estadd local space  					" "
		estadd local groupfe 					"Yes"
		estadd local countrypair 				"`countrypair'"
		estadd local maxdist 					"`xdist'"
		estadd local polynomial				"No"
		estadd local nmeandep				"`meandep'"
		local column_number = `column_number' + 1
		local column_name = "column" + "`column_number'"
	

	
*		xi: areg interdep  rdist1*  varofi  [aw=weight] ,a(groupfe) cluster(border)
		reghdfe interdep rdist1*  varofi [aw=weight], absorb(groupfe) vce(cluster groupfe country)
		eststo `column_name'
		estadd local space  					" "
		estadd local groupfe 					"Yes"
		estadd local countrypair 				"`countrypair'"
		estadd local maxdist 					"`xdist'"
		estadd local polynomial				"Linear"
		estadd local nmeandep				"`meandep'"
		local column_number = `column_number' + 1
		local column_name = "column" + "`column_number'"	
	
	
	
	}


}


label var interdep  " "		

#delimit ;
esttab column1 column2 column3 column4 column5 column6 column7 column8
using "`nametable'",  replace
label
b(3)
se(3)
star(* 0.10 ** 0.05 *** 0.01)
width(1.0\hsize)
nonotes
nolines
noconstant
mgroups(  "\textit{Dep. Var: Lowest risk crop is grown}" "\textit{Dep. Var: Highest risk crop is grown}"  , pattern(1 0 0 0 1 0 0 0 0)
prefix(\multicolumn{@span}{c}{) suffix(}) span erepeat(\cmidrule(lr){@span})) 
 keep(varofi )
order(varofi )
  stats(groupfe					space	polynomial		space	r2		N					countrypair				space	maxdist							space	nmeandep,
    fmt(0						0		0				0		3		0					0					0		0								0		0)
layout(@						@		@				@		@		@					@					@		@								@		@)	 	
labels("Border FE"				" "		"RD Polynomial"	" "		"R2"	"\# Observations"	"\# Country pairs"	" "		"Max distance to border (km)"	" " 	"Mean Dep. Var.")
 )	
noobs
booktabs
;
#delimit cr





